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We report new zero-field muon spin relaxation and neutron spin echo measurements in ferromagnetic 
(FM) (La,Ca)Mn03 which taken together suggest two spatially separated regions in close proximity 
possessing very different Mn-ion spin dynamics. One region corresponds to an extended cluster 
which displays 'critical slowing down' near Tc and an increasing volume fraction below Tc- The 
second region possesses more slowly fluctuating spins and a decreasing volume fraction below Tc- 
These data are discussed in terms of the growth of small polarons into overlapping regions of 
correlated spins below Tc, resulting in a microscopically inhomogeneous FM transition. 

PACS numbers: 75.30.Vn, 75.40. Gb, 76.75.+i., 25.40.Fq 
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After several years of study it is clear that the rich- 
ness in the temperature-composition phase diagram for 
Lai_2;Caa;Mn03 ||l[ is produced by a strong interplay be- 
tween the spin, charge and lattice degrees of freedom in 
these materials Q. Of particular interest has been the 
ferromagnetic (FM) insulator-to-metal transition and its 
accompanying large magnetoresistance. Millis and col- 
laborators p] first suggested that the theoretical descrip- 
tion of this transition must include the electron-phonon 
Jahn- Teller (JT) coupling, in addition to the double- 
exchange (DE) B spin-spin interaction, thus invoking 
polaronic degrees of freedom. Experimental signatures 
of small-polaron hopping have emerged from transport 
measurements above the FM critical temperature Tc [pi , 
leading to a theoretical description of magnetoelastic po- 
larons in terms of a small cluster of aligned Mn spins 
surrounding a local JT lattice distortion [g|J§]. 

Despite these important insights, an experimental 
and theoretical description of the FM transition encom- 
passing all three important degrees of freedom (lattice, 
charge, and spin) remains incomplete. The persistence 
of local lattice distortions below Tc is clear from neutron 
pair-distribution function (PDF) and x-ray absorption 
fine structure (XAFS) measurements 0,§|- A two-fluid 
model for the charge degrees of freedom is consistent with 
the evolution of localized charges into itinerant carriers 
near Tc P] . A comparable description of the behavior of 
the spin system near Tc is lacking, however. 

In this Letter we present new muon spin relaxation 
(//SR) and neutron spin echo (NSE) measurements in FM 
Lai_2:Caa;Mn03, x ~ 0.3. Neutron scattering has been 
used to study the low-temperature magnetic properties in 
FM perovskites flQ] , and has revealed a broad peak cen- 



tered around zero energy transfer which coexists with Mn 
spin waves near Tc |ll| ]. An unambiguous explanation 
for this peak is still lacking, however. The NSE technique 
|12| allows a direct measure of the spin-spin correlation 
function S{q,t) at much higher resolution and for longer 
correlation times than are achievable with conventional 
neutron scattering. The muon is a local probe (bound 
within about 1 A of an oxygen atom in oxide materials 
|l3[ ) and is, therefore, sensitive to spatial inhomogene- 
ity in spin fluctuation rates. Our /iSR and NSE results 
are consistent with a spatially distributed FM transition 
involving at least two interacting regions with very dif- 
ferent spin dynamics. We are able to follow the evolution 
of these different regions with /iSR over the temperature 
range 0.7 Tc <T < Tc- 

The /xSR data were taken at the Paul Scherrer Institute 
in Villigen, Switzerland, and at TRIUMF in Vancouver, 
Canada, between 10 and 300 K. The NSE data were 
taken at 275 K with the INll spectrometer at the Insti- 
tute Laue Langevin, using incident neutron wavelengths 
of 4 A and 6 A. All measurements were performed in 
zero applied field (< 1 Oe) on polycrystalline samples. 
(NSE measurements below Tc were impossible in zero 
field because the neutrons were completely depolarized 
by randomly oriented FM domains.) The /iSR sample 
[x = 0.33, Tc — 262(3) K from magnetization measure- 
ments] was from the same batch as that used in Ref. 
^. For the NSE sample x = 0.30 and Tc = 250(3) 
K. In general polycrystalline materials of (La,Ca)Mn03 
are more compositionally homogeneous than comparable 
volumes of single crystals, because of relatively less Ca 
evaporation and fluctuating growth conditions. 

The /xSR data could be fit to a relaxation func- 



tion Gz{t) = Gosc{t) + Grix(i), corresponding to os- 
cillating and relaxing terms, respectively. Goscit) is 
given by ^osc exp(-t/r2) cos(27ri/^< + 0^), where i^^ is 
the muon precession frequency, I/T2 is the inhomoge- 
neous damping rate, and 27ri'^ — j^B, with 7^ the 
muon gyromagnetic ratio and B the local magnetic field. 
Giix{t) was first fit to the stretched-exponential form 
A,-ixex.p[—{t/Ti)^], where 1/Ti is a characteristic spin- 
lattice relaxation rate and a value of .ftT < 1 implies a dis- 
tribution of rates. For our experiment Aqsc + ^rix ~ 0.2 
independent of temperature. 

The current /zSR data are of greater statistical pre- 
cision than previously reported [[l4| , and are taken at 
smaller temperature intervals near Tc, allowing a more 
refined interpretation of the relaxation phenomena. Fit 
parameters are shown in Figure 1. The observed zero- 
field frequencies i^^(T) tend to zero at a temperature of 
266 ±2 K, in agreement with Tq from the magnetiza- 
tion data mentioned above. The amplitudes Aqsc (not 
shown) indicate that the entire sample volume experi- 
ences growth of the sublattice magnetization. The higher 
statistical accuracy allows observation of the rapid de- 
cline of K to about 0.2 just below Tq, whereas previ- 
ously it was necessary to freeze K a,t a somewhat arbi- 
trary value of 0.5 below about 270 K |1J]. Because K 
and l/Ti are highly correlated, the l/Ti values reported 
here are different from those reported earlier. 

For rapid fluctuations the local muon relaxation rate A 
is given by A ex 7f,J2q\^B{q)\'^T{q), where \SB{q)\ is 
the amplitude of the fluctuating local field and T{q) is 
the Mn-ion correlation time. A distribution of A im- 
plies that SB and/or r are distributed but does not de- 
termine the distributions separately, whereas S{q, t) ob- 
tained from NSE measurements gives a direct measure 
of the distribution of T{q) only. Figure 2 displays S{q, t) 
ai T — 275 K in Lao.yoCao.aoMnOa. The data can be fit 
to S{q,t) = exp(— (i/r(q))'^), S{q,0) = 1; the parame- 
ters from these fits are given in Table I. These results, 
together with the /iSR data, directly confirm that the 
Mn-ion correlation times are spatially distributed. Note 
that /3 ~ 1 for g > 0.10 A"! and /3 < 1 for q < 0.10 k'^. 

The use of a stretched exponential to describe the 
/iSR and NSE relaxation functions is not unique and 
may not be physically appropriate, as we now dis- 
cuss. A stretched exponential relaxation function im- 
plies a continuous distribution P(A) of rates A defined 
by /dAP(A)exp(-Ai) = exp[-(At)"], which is broad 
for small a. Such a distribution is applicable to di- 
lute spin glasses, where a distribution of energy barri- 
ers yields a broad distribution of relaxation rates, and 
a given spin may relax via many possible channels p5[ . 
Thus NSE measurements in dilute spin glass CuMn |lfl] 
show no appreciable q dependence of S{q,t), an indica- 
tion of the many relaxation channels available to a given 
spin. This is in marked contrast to Lao.roCao.aoMnOa 
(Fig. 2), where the (/-dependence of S{q,t) is likely due 



to a spatial distribution of relaxation rates, as discussed 
below. We have therefore chosen to adopt a bimodal 
distribution P(A) ex AfS{X — A/) + As6{\ — As), corre- 
sponding to a two-exponential fit for S{q,t) and Giix(i). 
Further justification for this is given below. 

The results from fitting 

S{q,t) = AFiq)eM-t/rFiq)) + As{q)eyip{-t/'rs{q)), 
where Ap + As = 1, are given in Table I. The two com- 
ponents are labeled 'fast' (small r) and 'slow' (large r). 
The slow component sets in only for q < 0.10 A~^, 
and the correlation times Tp and ts differ by about 
an order of magnitude at the smallest q values. The 
/iSR data are also well described using a two-exponential 
form Grix(i) = ^/exp(— A/f) -I- Asexp(— A^t). Here 
Af + As ^ Aix- We have labeled the rates A/ and A^, in 
analogy to the NSE analysis. Again, the two exponential 
rates differ by at least a factor of ten, allowing good con- 
vergence of the fits. We note that non-uniqueness of the 
fitting function for both the NSE and fiSK data is a gen- 
eral property of sub-exponentially decaying curves [[l7[ . 
The justification for which function to use must therefore 
ultimately reside in a physical interpretation of the data. 

Figure 3 shows the temperature dependence of Af , Ag , 
A/, and As. For temperatures below 185 K the overall 
relaxation rate is too small to resolve two exponentials 
clearly, and above 273 K the exponent K approaches 1. 
Note that there is consistency between the NSE and /iSR 
measurements, taking into account that /iSR rates in- 
volve a sum over all q, which in ferromagnets is most 
heavily weighted for g ~ 0. The amplitudes of the fast 
and slow components near Tc are roughly equal. Al- 
though the ^SR data do not give precise values for Tf 
and Ts, because the SB(q) are not known, the fact that 
As/A/ and ts/tp (at small q) are both ^ 10 is quite 
consistent. 

The /iSR data in Fig. 3 exhibit the following important 
trends. The peak value of A/ coincides with v^{T) — > 
at r ~ Tc, unlike the stretched-exponential fits in Fig. 1, 
where the peak values of 1/Ti occur at distinctly lower 
temperatures (~ 255 K) than where v^{T) — > 0. Ordi- 
narily the spin-lattice relaxation rate is expected to peak 
at Tc, where the susceptibility and the spin-spin corre- 
lation times reach their maxima. This does not occur for 
the stretched-exponential fits because the rapidly chang- 
ing values of the exponent K near Tc result in a chang- 
ing functional form of Grix(r); then the derived values 
of 1/Ti are not consistent from one temperature to the 
next. Thus the two-exponential model function seems 
more appropriate. By contrast, the temperature depen- 
dence of As shows no significant maximum near Tc, but 
rises slowly below the critical temperature. Near Tc, Af 
is about 60% of the total relaxing amplitude and grad- 
ually increases as the temperature is reduced, reaching 
about 75% at T/Tc ^ 0.7. The dotted fines (Fig. 3, 
top) are guides to the eye, and indicate that the trend of 
Af is to approach 100% of the fluctuating amplitude at 



T = OK. 

The ^SR and NSE data suggest the following qualita- 
tive interpretation: there are two spatially separated re- 
gions in the sample, characterized in our measurements 
by very different Mn spin dynamics and temperature- 
dependent volumes. The muon relaxation samples these 
regions locally through its short-ranged dipolar cou- 
pling (matrix element oc r~^), predominantly to its two 
nearest- neighbor Mn spins. The /^SR rates involve sums 
over all q, but are dominated by FM fluctuations. NSE 
is sensitive to the spatial configuration of the spin fluc- 
tuations through the q dependence of the correlation 
times. The single-exponential relaxation observed for 
larger q reflects nearly single-ion dynamics which are 
not associated with the FM transition, whereas the two- 
component relaxation evident for smaller q shows that 
the FM dynamics are inhomogeneous. The length scale 
associated with this crossover is a few lattice spacings, 
which is roughly consistent with the 'droplets' found 
in small-angle neutron scattering in lightly-doped man- 
ganites iQ. We have also studied ^SR relaxation in 
undoped LaMnOa and CaMnOa, both of which ex- 
hibit temperature-independent single-exponential relax- 
ation above 140 K with rates 0.05-0.10 /is~^. This is 
very different than the behavior of either component 
in Lao.eTCao.asMnOa, which is strong evidence against 
large-scale separation of Mn'^+ and Mn'^^ ions. 

We postulate that the two observed relaxation com- 
ponents constitute the spin signature of magnetoelastic 
polarons as the system becomes FM. The more rapidly 
relaxing component, which shows a peak in the relax- 
ation rate at Tc (as expected for a FM phase transition), 
is attributed to spins inside the overlapping polaron re- 
gions. The volume fraction of these regions grows below 
Tc as the polarons continue to grow. The second compo- 
nent is characterized by much slower Mn-ion relaxation 
rates, shows no sign of critical slowing down, and occu- 
pies a diminishing volume fraction as the temperature is 
lowered. We associate this component with the spins at 
the boundaries of, or between, the polarons. 

In this picture the interior regions are characterized 
below Tc by relatively low (metallic) resistivity, an aver- 
age Mn valence between 3+ and 4+ (corresponding to the 
rapid motion of charge carriers via the DE interaction), 
and a smaller local JT lattice distortion than is found in 
LaMnOa. The rapid band- like carrier motion produces 
more rapid Mn spin fluctuations than in the less-metallic 
exterior regions, where spin and charge motion are frus- 
trated by more extreme local lattice distortions. These 
larger distortions limit the polaron size at a given tem- 
perature. Note that the two regions must be in close 
proximity to one another to explain the q dependence of 
S{q,t), and thus their spins should interact. This inter- 
action presumably causes weak FM alignment in the less- 
metallic regions, and may also explain the broad peak in 
A/ around Tc-, which is uncharacteristic of conventional 



critical dynamics. This breadth may also be caused by a 
distribution of polaron sizes (and local Tc values) . 

We compare these data to other measurements. Our 
picture is conceptually similar to the two-fluid model of 
resistivity Q , with one exception. The /zSR data indicate 
that slowly fluctuating spins in relatively less metallic re- 
gions persist far below Tc, whereas in the two-fluid model 
the growth of the conducting charge fraction reaches es- 
sentially 100% just below Tc- This can be explained 
by realizing that when a conducting path is reached (at 
Tc) the more resistive regions are short-circuited, even 
though they may still occupy a considerable volume frac- 
tion. Persistence of inhomogeneity below Tc is reinforced 
by local probes of the lattice structure ||3-i|j|9| 7 which are 
consistent with the gradual loss of structural inhomogene- 
ity below Tc, similar to that reflected in the spin- lattice 
relaxation rates reported here. Finally, it is likely that 
the central peak observed in inelastic neutron scattering 
|ll| corresponds to our slowly relaxing component. 

We have presented evidence for two time scales in the 
spin system of (La,Ca)Mn03, and associated these with 
a fine-scale spatially-distributed FM transition due to the 
formation and accretion of magnetoelastic polarons. Re- 
alistic theories of these slow inhomogeneous spin dynam- 
ics have yet to be developed. The spatial distribution 
of spin fluctuation rates is probably related to the dis- 
ordered spatial distribution of La and Ca ions and the 
corresponding local fluctuations of the local lattice dis- 
tortions, which influence the spin subsystem via magne- 
toelastic coupling. A theoretical model which includes 
disorder and couples JT pseudo-spin variables with the 
Mn-spin degrees of freedom may shed light on the mys- 
tery of the slow, inhomogeneous spin dynamics reported 
here. 
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FIG. 1. Temperature dependence of the muon precession 
frequency i^^ (top), average spin lattice relaxation rate l/Ti 
(middle) and exponent K (bottom) for the stretched expo- 
nential relaxation function Grix. 



FIG. 2. Spin-spin correlation function S{q, t) measured 
with NSE spectroscopy. The data for q > 0.10 A~ were 
taken at incident wavelengths of 4 A and 6 A and hence ex- 
tend to shorter times than data for q < 0.10 A^^, taken only 
at 6 A. 



FIG. 3. Temperature dependence of /iSR amplitudes 
Af/Arix and As/Aii^ (top) and their respective relaxation 
rates A/ (middle) and As (bottom) for a two-exponential fit. 
The subscripts / (filled symbols) and s (open symbols) refer 
to fast and slow Mn relaxation rates, respectively (see text). 
The dotted lines in the top frame show the trends extrapo- 
lated to r = K. Note the different temperature scales for 
the top and lower two frames. 



TABLE I. Fitting parameters for NSE spin-spin correlation 
function S{q,t) in Lao.roCao.aoMnOs, T — 275 K. Symbols 
are defined in the text. 
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